Spontaneous pattern formation in an acoustical resonator 
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A dynamical system of equations describing parametric sound generation (PSG) in a dispersive 
large aspect ratio resonator is derived. The model generalizes previously proposed descriptions of 
PSG by including diffraction effects, and is analogous to the model used in theoretical studies of 
optical parametric oscillation. 

A linear stability analysis of the solution below the threshold of subharmonic generation reveals 
the existence of a pattern forming instability, which is confirmed by numerical integration. The 
conditions of emergence of periodic patterns in transverse space are discussed in the acoustical 
context. 



I. INTRODUCTION. 



In recent years, pattern formation in systems driven far from equilibrium has became an active field of research in 
many areas of nonlinear science. Apart from peculiarities of particular systems, an outstanding property of pattern 
formation is its universality, evidenced when the dynamical models describing the different phenomena (either hydro- 
dynamical, chemical, optical or others) can be reduced, under several approximations, to the same order parameter 
equation. These equations are few and well known, such as Ginzburg-Landau or Swift Hohenberg [1]. 

A key concept in pattern forming systems is the aspect ratio. When the evolution of the variables is restricted 
to a bounded region of space, or cavity, the aspect ratio is defined as the ratio of transverse to longitudinal sizes 
of the cavity. In hydrodynamical Rayleigh-Benard convection, for example, the aspect ratio is determined by the 
ratio between the height of the fluid and the area of the container. In problems of nonlinear wave interaction, this 
parameter is related to the Fresnel number, usually defined as 



where a is the characteristic transverse size of the cavity (for example, the area of a plane radiator), A is the wavelength 
and L is the length of the cavity in the direction of propagation, considered the longitudinal axis of the cavity. 
Spontaneous pattern formation is observed in large aspect ratio nonlinear systems driven by an external input, 
, where the possibility of excitation of many transverse modes (a continuum for an infinite transverse dimension) is 
considered. When the amplitude of the external input reaches a critical threshold value, large enough to overcome the 
losses produced by dissipative processes in the system, a symmetry breaking transition occurs, carrying the system 
from an initially homogeneous to a inhomogencous state, usually with spatial periodicity. These solutions have been 
often called dissipative structures [2]. A paradigmatic example in pattern formation studies has been the Rayleigh- 
^ ■ Benard convection in a fluid layer heated from below, where roll or hexagonal patterns are excited above a given 
temperature threshold [1]. 

This scenario differs with that observed in small aspect ratio systems such as, for example, a waveguide resonator 
of finite cross section. In this case, boundary-induced spatial patterns arc selected not by the nonlinear properties of 
the cavity, but by the transverse boundary conditions, and correspond the excitation of one or few transverse modes 
of the cavity [3]. In this sense, boundary-induced patterns are of linear nature. 

Guided by the analogies with other physical systems, we can expect that an ideal system for such effects to be 
observed in acoustics consists in a resonator of plane walls (acoustical interferometer) with infinite transverse size. In 
practice, the large aspect ratio condition could be fulfilled if the pumped area is finite, but large in comparison with 
the spatial scale imposed by the cavity length and the field wavelength, as follows from (1). 

Parametrically driven systems offer many examples of spontaneous pattern formation. For example, parametric 
excitation of surface waves by a vertical excitation (Faraday instability) [6], vibrated granular layers [7], spin waves 
in ferrites and ferromagnets and Langmuir waves in plasmas parametrically driven by a microwave fields [8], or the 
optical parametric oscillator [9,10] have been studied. 

The behaviour of nonlinear waves in large aspect ratio resonators has been extensively studied in nonlinear optics 
(for a recent review, see [4]), where a rich variety of patterns has been observed. On the other side, it is well known 
that optical and acoustical waves share many common phenomena, under restricted validity conditions [5]. 
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In nonlinear acoustics, a phenomenon belonging to the class of the previous examples is the parametric sound 
amplification. It consists in the resonant interaction of a triad of waves with frequencies u>o, ui\ and u>2, for which the 
following energy and momentum conservation conditions arc fulfilled: 

CJ = Wi + 0J2, 

k Q = fcj + k 2 + Ak, (2) 

where Ak is a small phase mismatch. The process is initiated by an input pumping wave of frequency u> which, 
due to the propagation in the nonlinear medium, generates a pair of waves with frequencies lo\ and ui2- When the 
wave interaction occurs in a resonator, a threshold value for the input amplitude is required, and the process is called 
parametric sound generation. In acoustics, this process has been described before by several authors under different 
conditions, either theoretical and experimentally. In [11-13], the one dimensional case (colinearly propagating waves) 
is considered. In [3], the problem of interaction between concrete resonator modes, with a given transverse structure, 
is studied. In both cases, small aspect ratio resonators containing liquid and gas respectively are considered. More 
recently, parametric interaction in a large aspect ratio resonator filled with superfluid He 4 has been investigated 
[16,17]. 

The phenomenon of parametric sound generation is analogous to optical parametric oscillation in nonlinear optics. 
However, an important difference between acoustics and optics is the absence of dispersion in the former. Dispersion, 
which makes the phase velocity of the waves to be dependent on its frequency, allows that only few waves, those 
satisfying given synchronism conditions, participate in the process. 

In a nondispcrsivc medium, all the harmonics of each initial monochromatic wave propagate synchronously. As 
a consequence, the spectrum broadens during propagation and the energy is continuously pumped into the higher 
harmonics, which eventually leads to shock formation. 

Optical media are in general dispersive, but acoustical media not. In finite geometries, such as waveguides [15] or 
resonators [14], the dispersion is introduced by the boundaries. Different cavity modes propagate at different angles, 
and then with different "effective" phase velocities. However, in unbounded systems boundary- induced dispersion is 
not present. 

Other dispersion mechanisms have been proposed in nonlinear acoustics, such as bubbly media or layered (periodic) 
media [18]. In all this systems, dispersion appears due to the introduction of additional spatial or temporal scales in 
the system, which makes sound velocity propagation to be wavelength dependent. Other proposed mechanisms are, 
for example, media with selective absorption, in which selected spectral components experience strong losses and may 
removed from the wave field [19]. 

Pattern formation in acoustics has been reported previously in the context of acoustic cavitation, [24,25] , where the 
coupling between the sound field amplitude and the bubble distribution is considered. We do not consider here the 
dynamics of the medium, which is assumed to be at rest, but the coupling between different frequency components. 

The aim of the paper is twofold. On one side, a rigurous derivation of the dynamical model describing parametric 
interaction of acoustic waves in a large aspect ratio cavity is presented. The derived model is analogous to the system 
of equations describing parametric oscillation in an optical resonator, and consequently their solutions are known. On 
the other side, among other peculiarities, the model predicts a pattern forming instability, which is confirmed by a 
numerical integration. We review these properties in the acoustical context, giving some estimations of the acoustical 
parameters which can make the model closer to real operating conditions. 



II. MODEL EQUATIONS. 

The starting point of the theoretical analysis is the nonlinear wave equation, which written in terms of the pressure 
p reads 

2 1 d 2 P 8 d 3 P _ _ e d 2 p 2 _ ( 2 1 a 2 \ 

where C = \pqv 2 — 2 V PQ ^ ^ s * ne lagrangian density, with po the ambient pressure and v the particle velocity, 5 is the 
diffusivity of sound, defined as 

*-s(i' +<+ "(s-i))' (4) 

and e = 1 + B/2A is the nonlinearity coefficient. 
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In the case of colinearly propagating plane waves, the lagrangian density term vanishes, since the linear impedance 
relation p = pocv holds. However, even in the case of slightly diverging waves (plane waves propagating at a small an- 
gle), the last term in (3) is much smaller in magnitude than the other terms in the equation, and in this approximation 
its effects can be ignored [15]. 

Furthermore, in this case the nonlinearity parameter can be considered to be independent of the interaction angle. 
It has been shown in [15] that, for the process described by Eqs.(2), when the waves with frequencies uu\ and cj 2 
propagate at an angle 9, the nonlinearity parameter for the u)q wave is expressed as 



which reduces to e = 1 + B/2A when 9 is small. 

Under this assumptions, the field distribution can be properly described by the wave equation 

2 _ \_cPp_ _ 5_(Pp 

p c 2 dt 2 c 4 at 3 p c 4 at 2 ' 1 ' 

which is the well known Westervelt equation. 

Equation (6) describes the propagation of waves in a nondispersivc homogeneous medium, but also in a medium 
that possess some dispersion mechanism, such as a bubbly liquid when the field frequencies are much lower than the 
resonance frequency of the bubbles [20]. It must be noted that, in a dispersive medium, the nonlinearity parameter e 
depends on the dispersion mechanism. 

We consider in the following that the wave interaction is only effective among three resonant frequencies, for 
which the relations (2) hold. Taking into account that, due to reflections in the walls, there exist waves propagating 
simultaneously in opposite directions, the field inside the resonator can then be expanded as 

3 

p(r,i) = ^P,(r,t), (7) 

where Pj(r,t) are the wave components related with the frequency Uj, given by 

Pj{r,t) = pj(x,y,t) cos (hp) e"^* + c.c. (8) 

where fc| = mir/L is a cavity eigenmode. 

A more general solution can be proposed, consisting in a superposition of quasi-planar waves, the field being 
decomposed in forward (F) and backward (B) waves, respectively. The quasi-planar assumption implies that the 
amplitudes may depend on the longitudinal coordinate z, representing a slow evolution compared with the scale given 
by kj. The field at frequency uij is then expressed as 

Pj( r ,t) = [Fj(r,t)e ik i z + B J (r,t)e- lk > z ] e"^* + c.c. (9) 

where r = (x, y, z) . 

The case w = 2uj\ corresponds to degenerate interaction (lo\ = u 2 = ojq/2), and describes the process of second 
harmonic generation or subharmonic parametric generation, depending on whether the pumping wave oscillates at ujo 
or u)\. In the following, the degenerate interaction case is considered, where 

p(r,t) = P (r,t) + P 1 (r,t) (10) 

with the amplitudes given in (9). 

Substitution of (10) in (6), and projecting the resulting equation on each of the mode frequencies, two coupled wave 
equations are found, 

V 2 P 1 d*P 3 5j d*Pj e d 2 (P 2 ) 3 

3 c 2 dt 2 + c] at 3 PQ c) at 2 [ ' 

where j = 0, 1 and Cj — uoj/kj = c (u)j) is the phase velocity of the waves. 

The evolution equations for the amplitudes pj at the different frequencies can derived be from (11) under several 
approximations. The complete derivation is given in the Appendix. If we assume that (i) the amplitudes are slowly 
varying in z and t, (ii) the reflectivity IZj at the boundaries is high, (Hi) the field frequencies ojj are close to one 
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resonator frequency, uij and that (iv) only one longitudinal mode is excited by each frequency component, the evolution 
is ruled by the following dynamical equations: 

=£;-7o(l + «Ao)po + «aoV_ L po-^Pi, 

= -7i (1 + i A i) Pi + iai Vjj>i - i*jP*iPo (12) 

where (x,y,t) are the amplitudes corresponding to (8), the nonlinearity parameter is denned as = ekj/2poCj, 
and the asterisk denoted complex conjugation. The other parameters in (12) are the pump E, detuning Aj, losses jj 
and diffraction aj and are defined as 

E=^^Ye^, (13) 




7i = l^ + *. (15) 

* ^ I" < 16 > 

where F is the amplitude of the incident wave at frequency uj , TZj arc the reflectivities at the boundaries, 9j = 
2L (ujj — ujj) /cj is the frequency mismatch, and Kj — Skj/2 is a loss factor related with the diffusivity of sound. 
Note that the total loss parameter 7^ takes into account both loss mechanisms, due to absorption and reflectivity 
at the boundaries (energy leakage of the resonator). Also, from (14) and (15) the detuning can be written as 

Finally, the field amplitudes are normalized to leave the model in its final form, 



1 dA 



= £-{!+ i A )A Q + ia V 2 A - A\ 



70 dt 
1 dA 



71 dt 



1 - = -{l + iA 1 )A l +ia l V 2 A l +AuA{. (17) 



where Aq = i(bi/2)p , A\ = i-^/bobi/Spi and £ = i(bi/2)E. 

The system of equations (17), together with their complex conjugated, consists in the model for parametric interac- 
tion of acoustic waves in large aspect ratio resonators. Eqs.(17) are suitable for the description of the spatio-temporal 
evolution of the pressure envelope waves oscillating at the fundamental (w ) and subharmonic (u)\) frequencies. 

These equations have been studied in the context of nonlinear optics, as a model for optical parametric oscillation 
[9]. In the following we review some of the basic properties of the solutions of Eqs.(17), and their application to the 
case of the acoustic resonator. 

III. HOMOGENEOUS SOLUTIONS 

The stationary and spatially homogeneous solutions of (17) are obtained when the temporal derivatives and the 
transverse diffraction term vanish. In this case, the model reduces to the one derived for parametric sound generation 
in a resonator with rectangular cross section a x b, described in [3], where its homogeneous solutions were obtained. 
We next review these solutions in the present notation. 

The simplest homogeneous solution corresponds to the trivial solution, 

1 + «A 

characterized by a null value of the subharmonic field inside the resonator. This solution exists for low values of the 
pump amplitude (below the instability threshold to be discussed in the next section). 
For larger pump values, also the subharmonic field has a nonzero amplitude, given by 



I A 1 1 2 = -1 + A Ai ±\Js 2 — (Ai + AO 2 , (19) 
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while the stationary fundamental intensity |^4 | takes the value 

|A | 2 = 1 + A?, (20) 

which is independent of the value of the injected pump. 

The emergence of this finite amplitude solution corresponds to the process of parametric generation. In the frame of 
the plane wave model it has been shown theoretically and confirmed experimentally [3], that the trivial solution (18) 
bifurcates in the nontrivial one, and the subharmonic field emerges with an amplitude given by (19). This bifurcation 
is supercritical when AoAi < 1, and subcritical when AoAi > 1. In the latter case, both homogeneous solutions can 
coexist for given sets of the parameters. 

In the next section this scenario is generalized by including diffraction effects in the model. 



IV. LINEAR STABILITY ANALYSIS 



In order to study the stability of the trivial solution (18) against space-dependent perturbations, consider a deviation 
of this state, given by 



Aj (x, y, t) = Aj + SAj (x, y, t) . 



(21) 



Assuming the deviation to be small, after substitution of (21) in (17) the resulting system can be linearized in the 
perturbations 5Aj . This leads to the linear system of equations 



1 dSA 
7o dt 
1 88A l 



(1 + iA )SA + ia V 2 SA , 
(1 + iAi)<L4i + ia 1 V 2 5A 1 - AqSA*. 



7i dt 

The generic solutions of (22) are of the form 

(SA^SA*) oc e *( fc -0* e ik 



(22) 



(23) 



where A(kj_) represents the growth rate of the perturbations, and is the transverse component of the wavevector, 
which in a two-dimensional geometry obeys the relation |kj_| 2 = k 2 + k 2 . 

Substitution of (23) in (22) and its complex conjugates, written in matrix form, allows to evaluate the growth rates 
A as the eingenvalues of the stability matrix L. This is given by 



T iu \ \ 







(24) 



where the block matrices arc defined as 



V o o Ll {k±) J 

, . / -1 + iA + ia k\ 
MM=7o^ o -l- 2 Ao-za £d 



and 



Li (k±) = 7i 



-1 + iAi +iaik\ y^f^ 

, , £ A -1 — iAi — iaik 2 . 



(25) 



(26) 



The eigenvalues of Lq have always a negative real part, as follows from (25), and do not predict any instability. 
The eigenvalues of L\ are associated with the instability which gives rise to the subharmonic field, and are given by 



A± (fc ± ) = -l± 



£ 2 



2 \ 2 



1 + Al 



- {A 1 + a 1 k 2 ± ) 



(27) 



Note that only the eigenvalue with the positive sign A+ can take positive values, and reach the instability condition 
Re (A) > 0. 
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The eigenvalue is wavenumber dependent, which means that not all the perturbations in the form of transverse 
plane-wave modes (23) grow at the same rate. The maximum growth rate follows from the condition d\+/dk = 0. 
Two different cases, depending on the sign of the subharmonic detuning, must be considered: 

If Ai > 0, which corresponds to a subharmonic frequency larger than that of the closest cavity mode, the eigenvalue 
shows a maximum at 

k± = 0. (28) 

In this case, the emmited subharmonic wave travels parallel to the cavity axis, without spatial variations on the 
transverse plane. The solution remains homogeneous, and its amplitude is given by (19). 

On the contrary, if Ai < 0, corresponding to a field frequency smaller than that of the cavity mode, the maximum 
of the eigenvalue occurs at 



The corresponding solution is of the form (23), which represents a plane wave tilted with respect to the cavity axis. 
This solution presents spatial variations in the transverse plane, and consequently pattern formation is expected to 
occur. The two kind of instabilities are represented in Fig.l, where the eigenvalue as a function of the wavenumber is 
plotted, in the cases of positive (a) and negative (b) detuning, for given values of the pump above threshold. 

The value of the transverse wavenumber given by (29) can be interpreted in simple geometrical terms: kj_ corre- 
sponds to the tilt of the wave necessary to fit the longitudinal resonance condition. This is a linear effect related with 
diffraction, and it is represented in the figure 2. 

This can be also analytically shown by inspecting the solutions of the dyn amical eq uations for a transverse wave 
in the form (23): the frequency of the cavity mode is given by u> = c|k| = cy/k\ + fcf. Since, in the small detuning 
case, the relation k± « k z holds, it can be written approximately by 

ck z + ^-ki=ck z + Aoj (30) 
2k z 

where Alo is the transverse contribution to the mode frequency This contribution arises from the diffraction term 
ia\7 2 A — i^-V 2 A in the dynamical equations. 

Since k± is the modulus of the wavevector, the linear stability analysis in two dimensions predicts that a continuum 
of modes within a circular annulus (centered on a critical circle at |kj_| = k± in (k x , k y ) space) grows simultaneously 
as the pump increases above a critical value. This double infinite degeneracy of spatial modes (degenerate along a 
radial line from the origin and orientational degeneracy) allows, in principle, arbitrary structures in two dimensions. 

The threshold (the pump value at which the instability emerges) depends also on the sign of the detuning, and is 
obtained from the condition Re A = 0. From (27) it follows that 




£ th = ^/l + A§^/l + (A 1 +a 1 fc 2 ) 2 . (31) 

From the previous analysis, again two cases must be distinguished. For positive detunings, the mode homogeneous 
k± = experience the maximum growth, and the threshold occurs at a pump value 

Sth = /l + Ag^l + A?, (32) 

which is the same found in [14]. The solution above the threshold is given by (19). 

For negative detunings, the instability leads to a non homogeneous distribution, with a characteristic scale given 
by the condition a\k\ = — Ai. The corresponding threshold for such modes is 



£ th = yJl + *$. (33) 



Besides the existence of a pattern forming instability, a relevant conclusion of the previous analysis is the prediction 
of a decrease in the threshold pump value of subharmonic generation when diffraction effects are included [compare 
Eqs.(32) and (33)]. This fact is specially important in acoustical systems where the nonlinearity is weak, since in this 
case the instability threshold, appearing for high values of the driving, can be lowered by a factor of yl + Aj. 

The predictions of the stability analysis correspond to the linear stage of the evolution, where the subharmonic field 
amplitude is small enough to be considered a perturbation of the trivial state. The analytical study of the further 
evolution would require a nonlinear stability analysis, not given here. Instead, in the next section we perform the 
numerical integration of Eqs.(17), where predictions of the acoustic subharmonic field in the linear and nonlinear 
regime are given. 
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V. NUMERICAL SIMULATIONS 



In order to check the analytical predictions of the linear stability analysis, we integrated numerically the system (17) 
by using the split-step technique on a spatial grid of dimensions 64x64. The local terms, either linear (pump, losses 
and detuning) and nonlinear, are calculated in the space domain, while nonlocal terms (diffractions) are evaluated in 
the spatial wavevector (spectral) domain. A Fast Fourier Transform (FFT) is used to shift from spatial to spectral 
domains in every time step. Periodic boundary conditions are used. 

As a initial contition, a noisy spatial distribution is considered, and the parameters are such that a pattern forming 
instability is predicted (Fig. 3a). 

For small evolution times, the amplitude of the subharmonic remains small, corresponding to the the linear stage of 
the evolution. As follows from the linear stability analysis, an instability ring in the far field (transverse wavenumbcr 
space) is observed (Fig. 3b). 

For larger evolution times, the nonlinearity comes into play. In the nonlinear stage of evolution, a competition 
between transverse modes begins, and pattern selection is observed. In Fig. 4. it is shown a transient stage, where a 
labyrinthic pattern is formed. The final state, not shown in the figure, corresponds to the selection of a discrete set 
of transverse wavevectors, asymptotically resulting in a periodic pattern. 

VI. ACOUSTICAL ESTIMATES 

The theory presented in this paper has been suggested by analogies with nonlinear optical resonators, where good 
agreement between theory and experiment has been shown [4]. Next we estimate the required physical conditions to 
make the theory applicable to an acoustical nonlinear resonator. 

The mean field assumption implies that the field envelopes changes little during a roundtrip of the wave in the 
resonator. Acoustic waves, specially in the ultrasound regime, experience strong losses during propagation in the 
medium. It is then required the resonator to be short enough in order to avoid strong absorption. 

The longitudinal size of the resonator imposes a condition on the transverse scale, in order to maintain the Fresnel 
number large. For example, for a pump wave of 1 MHz in water the corresponding wavelength is A w 1500/10 6 s=s 1.5 
mm. If the resonator length is 2 cm and the walls (transverse size of the transducer) are squared with 10 cm each 
side, the Fresnel number gives F = 167, which can be considered in the large aspect ratio limit. 

With this geometry, the experimental setup proposed in the first observation of parametric sound generation in a 
liquid filled resonator [11] seems a good candidate for the observation of the predicted phenomena. In particular, it was 
shown in [11] that, close to the threshold of parametric generation, only the fundamental and subharmonic frequencies 
where present inside the resonator, in agreement with our theory. The resonator was formed by a PZT transducer, 
driven at a pumping frequency w Q , and a reflector placed at a variable distance L. The threshold was achieved when 
the transducer was driven at a voltage of nearly 300 V. The discrete spectrum, containing the fundamental w and the 
subharmonic co/2 frequencies, existed only close to the threshold (the relative pumping level being ldB). In some cases, 
pairs of signal-idler waves (corresponding to the nondegenerate process) were excited. For larger driving voltages, 
and even larger number of cavity modes, at different combination frequencies, appeared due to the recurrent process 
of parametric amplification. The evolution of the fields at such high pump values can not be predicted by the model 
proposed in this paper, which was derived under the assumption that the energy exchange occurs only between two 
modes. Although the spectrum is still discrete in this case, the large number of modes makes the spectral approach 
inconvenient for the theoretical study. Instead, the field approach, typically used in problems of nonlinear acoustic in 
nondispersive media, seems more convenient. A discussion concerning the applicability of the two approaches to the 
description of nonlinear wave processes in acoustics can be found in [26] . 

VII. CONCLUSIONS 

In this paper, the problem of spontaneous emergence of patterns in an acoustical interferometer is studied from the 
theoretical point of view. A model for parametric sound generation in a large aspect-ratio cavity is derived, taking 
into account the effects of diffraction. It is shown that the subharmonic field can be excited, when a threshold pump 
value is reached, characterized by a non-uniform distribution in the transverse plane, or pattern. This occurs when 
the field frequency is tuned below the frequency of the closest resonator mode (negative detuning). On the contrary, 
an on-axis field, with homogeneous distribution, is emmited. Traditionally, patterns arise from the imposition of 
external constraints (waveguiding) . The field then oscillates in modes of the resonator. The spontaneous emergence 
of patterns considered in this article presupposes no external transverse mode selection mechanism, but instead allows 
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the system to choose a pattern through the nonlinear interaction of a tipically infinite set of degenerate modes. The 
pattern formation process described here is related with the competing effects of nonlinearity and diffraction, and 
presents many analogies with similar systems studied in nonlinear optics, such as the two-level laser or the optical 
parametric oscillator, to which the model derived in this paper is isomorphous. 
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APPENDIX 

In this appendix we derive the model for parametric sound generation given in (17). The approach is based on the 
method proposed in [22] to describe the dynamics of laser fields in Fabry- Perot resonators. 

The first assumption is that the longitudinal variations of the fields are mainly accounted for by the plane wave, in 
which case the envelope amplitudes Fj and Bj can be considered as slowly varying functions in z and t. This leads 
to the inequalities 



d 2 F, 



dz 2 



« 



d 2 Fi 



dt 2 



« 



0F~ 



(A-l) 



which allow to neglect the second order derivatives in z and t. 

In the same way, applying the second relation of (A-l) to the dissipation term, it reduces to 



c] dt 3 



■&3 3 



(A-2) 



The condition (A-l) leads to the well known parabolic or eikonal approximation, in which the d'Alembertian 
operator acting on the pressure on the left-hand side of (11) can be approximated by 



cjdi 2 



( d 


1 d\ 


( d 2 


d 2 \ 


Ur 


C] dt) 


\dx 2 


V dy 2 ) 



(A-3) 



With (A-2) and (A-3) the wave evolution can be written as 



dt + 2iL v - lFj 



dBi dBi 



dz 



dt 



+ ^-VlB j -K j B j )=-ib j [p 2 ]', 



(A-4) 



where «;,• 



8-k 2 

is a parameter that accounts for dissipation, bj 



ekj/2poCj is a nonlinearity parameter, and 



where [p 2 ]j contains the terms in p 2 that oscillate with frequency — u)j, and takes into account slow and fast spatial 
variations. 

The nonlinear terms at the frequencies oj and u)\ are, respectively, 



-B 2 e- tkoZ + -B-tFx + -F 2 e lkoZ 
4 2 4 



\ B * Bu< 



— ik\ z 



1 p* D —'Aik-iz 



r>* tp Jiikiz 



-F*F e iklZ 
2 1 



(A-5) 



In order to eliminate the explicit dependence of the exponential factors in (??), a projection on two longitudinal 
modes is performed, multiplying by exp (±ikjz), and integrating over a full wavelength. This leads to a system of 
equations where the fields Fj and Bj are decoupled in the linear part (l.h.s), i.e., we get 



dF dF Cj 2 bj 1 r r , r 



<d{kjZ) 



(A-6) 
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for the forward waves, and 



for the backward waves. 

Substituting the nonlinear terms (A-5), only one of the contributions survive, leading to a couple of equations for 
each frequency, 

dF dF c 2 6 F 2 



for the fundamental wave, and 



- C0 ^7 ~dT 2^ v ^° - KaBa = (A " 8) 

-^ + ^ + ^^-^-4^ (A-9) 

for the subharmonic. 

Besides the existence of the counter-propagating waves, the cavity also imposes a condition which relates the 
amplitudes of forward and backward in given points of the resonator. The fields obey the following boundary conditions 
in the cavity: in z = 

F (x, y, 0, t) = VT^Bo (x, y, 0, t) + %Y, 

F 1 {x, y, 0, t) = y/n x B x (x, y, 0, t) , (A-10) 

and in z = L 

B a {x, y, L, t) = ^e- i6 °F (x, y, L, t) , 

Bi(x,y,L,t) = \fR\e~ idl F\ (x, y, L, t) (A-ll) 

where L is the length of the cavity, T$ is the transmitivity of the boundary to the input wave (pump), with amplitude 
given by Y, IZj is the reflectivity of the field pj at the boundary. Finally, 9j are the detunings (frequency mismatch 
with respect to the cavity), given by 

6o = 2wm - 2k L = ^ ~ U T ° , (A-12) 

6>i = 2nn - 2kyL = ^ ~ Ul , (A-13) 

being w| the frequency of the cavity (eigenmode) nearest to the field frequency ojj . 

The field evolution is described by the system of equations (A-8) and (A-9) together with (A-10) and (A-ll). 
However, the description can be greatly simplified after the introduction of the following changes [22]: 

Fo = e £(z-L)(lnKo-ie ) Fo + ^ e ^J_ {z _ L)y 

V 7?-o 2L 

Bo = e-^^-^Boe^ - ^L** -L (z — L)Y 

\/1Zq 2L 

F 1 = e TL^-m^i-i9l)p^ ( A _ 14 ) 

B x = e-^ iinni - Wl) B 1 e 1 ^ , 
The great advantaje of the latter changes is that the boundary conditions for these new fields are 

F o (0,t) =B (0,t), 

h{L,t) =Bi(L,t) (A-15) 
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which correspond to those of an ideal cavity with perfectly reflecting boundaries. This fact will be used later for 
determining the unknow longitudinal distribution. 

Substituting the new amplitudes in the evolutions equations we find the new system 



dF dF \\n1Z \+i6 ~ \lnKo\+i0 V% v t ?° , 

-dT +c ^ = - c ° — 2L — Fo+c ° & {z - L) vm Ye 2 + 

CO V%„.i e -£ * rff „p Mr.,.™ 



-- Yel * - 7T^ V ^ - k q F - z^V (z) Ft (A-16) 



2iV^o 2ifc 4 

where we have defined the exponential term 

V (z) = e 2r( z - L )( ln ^o-J0 o ) e ^(z-L)(ln 7^-10!)^ (A- 17) 

Consider now two additional conditions: first, that the reflectivity at the boundaries is high (weal losses), which is 
expressed as IZq — ► 1, and consequently 7$ — > 0. Second, that the field frequencies are close to a resonant frequency 
of the cavity, so that 8a — > 0. Under these conditions, also known as the mean-field limit in nonlinear optics, the 
exponential factor T> (z) approaches unity, and the second term at the l.h.s. of (A-16) can be neglected. Furthermore, 
in this limit the field amplitudes Fj and Bj approach to their real values, Fj and Bj. The equation (A-16) takes then 
a simplified form, which can be more conveniently written as 

ifr + = E - 70 (i + zAo) F ° + ia ° v2F ° - (A " 18) 

where we have defined the new pump E, detuning A , losses 70 and diffraction ao parameters as, 

e =?S/'- <A - i9) 

I In I c,- , , „„. 

7j= 2L +Kj ' (A " 21) 

a ' = ^ (A " 22) 

where j = 0, 1. Defined in this way, the detuning is a quantity of the order of 1. Also, 70 represents a measure of the 
lincwidth of the cavity modes. 

We have performed the derivation for the forward component of the fundamental wave. A similar procedure leads 
to the evolution equations for the other components, 

lit + c °lt = E ~ 70 (1 + lAo) B ° + maV2Bo ~ l i B ^ ( A " 23 ) 

fii? 3F b 

-df+ c ^ = -^( 1 + zA i) F i + ^iV 2 Fx - ijF*F , (A-24) 
+ Ci -q± = -71 (1 + iAO B 1 + ia l V 2 B l - i^B{B (A-25) 

The equations still keep a explicit z dependence. Consider now the new boundary conditions, given by (A-15). 
These conditions differ from that of the original fields in that they represent an ideal (lossless) cavity, which allows 
to express the field inside the cavity by means of the Fourier expansions 



00 

F i ( r > *) = XI p f ] ( x ' y ' *) cxp 0" 



rnrz 



n— — 00 



J' 



4- (r, t) = £ p$ n) (x, y, t) exp (-^) (A-26) 
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From (9), it follows that the total field with frequency ujj can be also written as 



oo 



Pj(r,t)= J2 2^" ) (x, 2/ ,t)cos[( 



ki + ™)z]e 



(A-27) 



n— — oo 



If we finally assume that the intermode spacing c/2L is large, we can consider that only the n = mode can be 
excited, and Fq and Bq became spatially uniform along the longitudinal axis of the cavity. In this case, the field can 
be described by (8), which is the usual description for waves in acoustic resonators. 

This last assumptions makes the dynamical evolution to be independent on z, and Eq. (A-18) take the final form 



which are the equations (12) given in the text. 
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dpo 
dt 

dpi 
dt 



E - 7o (1 + iA ) p + ia X7 2 p - i-^pj 
-71 (1 + iAi)pi + iaiV 2 pi - i^p\po 



(A-28) 
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Figure captions: 



Figure 1: Eigenvalue as a function of the perturbation wavenumber, for Ai = 1, E = 2.25 (a), and for Ai = — 1, 
E= 1.35(b). 



Figure 2: Schematic representation of the longitudinal cavity resonance condition for negative dctunings. 



Figure 3: The linear stage of the evolution. It is shown the pressure amplitude (left column) and the corresponding 
spectrum (right column) in the transverse plane. Pictures were taken at times t = 1 (a) and t = 5 (b) The parameters 
used in the integration are Ai = —2, Aq = — 1, E = 2, 70 = 71 = 1, a\ = 0.001 and ao = 0.0005. 



Figure 4: Nonlinear evolution and pattern selection. The same parameters as in Fig. 3 have been used. Pictures 
were taken at times t = 25 (a), t = 750 (b) and t — 1000 (c). 
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